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The scattering amplitude from a set of discrete states coupled to a continuum became known as 
the Fano profile, characteristic for its asymmetric lineshape and originally investigated in the context 
of photoionization. The generality of the model, and the proliferation of engineered nanostructures 
with confined states gives immense success to the Fano lineshape, which is invoked whenever an 
asymmetric lineshape is encountered. However, many of these systems do not conform to the initial 
model worked out by Fano in that i) they are subject to dissipative processes and ii) the observables 
are not entirely analogous to the ones measured in the original photoionization experiments. In this 
letter, we work out the full optical response of a Fano model with dissipation. We find that the 
exact result for the excited population, Raman, Rayleigh and fluorescence emission is a modified 
Fano profile where the typical lineshape has an additional Lorentzian contribution. Expressions to 
extract model parameters from a set of relevant observables are given. 


In a set of seminal papers spanning from 1935 to 1961, 
Beutler Q], Fano m and Friederichs [3] laid the basis of 
the theory to describe the absorption lineshapes of atomic 
photoionization experiments. These lineshapes present 
marked asymmetries which could not be explained by a 
simple Lorentzian resonance. The explanation was at¬ 
tributed to an interference between two photo-ionizaton 
pathways: one where the atom is ionized directly from its 
ground state and one where it is first excited to a higher 
discrete state which then ionizes (auto-ionizing states). 
The minimal Fano model consists in a discrete excited 
state coupled to a continuum set of excited states, both 
types of states being reachable by photo-excitation from 
the ground state. The resulting photo-fragmentation 
cross-section as a function of the excitation laser fre¬ 
quency lol is known as the Beutler-Fano or Fano profile: 


where q is the ratio of the transition dipole moment of the 
ground-discrete and ground-continuous transitions, and 
e = (Hujl — E e )/7 where E e is the energy of the discrete 
state relative to the ground state and 7 = mrV 2 is the 
linewidth of the excited state, induced by its coupling 
(per unit of energy) nV 2 to the continuum set of states, 
n being the density of states. 

Since the original photoionization experiments, the 
Beutler-Fano profile has been observed in an ever in¬ 
creasing variety of physical systems, and in particular 
in nanoscale structures mm- These include plasmonic 
nanostructures EH3, quantum dots, decorated nanopar¬ 
ticles IO! and spin filters tm. to name a few. Although 
the Fano theory was built on a scattering framework 
where the observable was the population on the contin¬ 


uum (i.e. the ionized electrons), the result continued to 
be applied (with remarkable success) to dissipative, non¬ 
scattering systems where the observable was not always 
the population in the continuum of states. As noted by 
A. E. Miroshnichenko in 2010, “a suitable theory for a 
quantitative description of these cases is still lacking” [5]. 

Some attempts to supplement the Fano model with 
dissipation processes have already been made. First by 
Fano himself in 1963 [T5] in a less-known paper, where 
the objective was to model pressure lineshape broaden¬ 
ing by atomic collisions. In this paper, only pure dephas¬ 
ing was considered which amounts to considering solely 
elastic collisions. In the eighties, Agarwal et al. [13] in¬ 
cluded population relaxation from the discrete excited 
state to an additional discrete state to model the compe¬ 
tition between atomic photo-ionization and spontaneous 
emission. K. Rzyzewski and J. H. Eberly |14j considered 
pure dephasing using a Wiener-Levy stochastic process 
to describe phase fluctuations. More recently Kroner et 
al. [T5] and Zhang et al. [IB] included a general dissipation 
process for a system of semiconductor quantum dots but 
the equations were solved approximatively by neglecting 
the population on the continuum set of states. It is worth 
noting that these three last works focus on intense field 
effects, which will not be adressed in the present letter 
and will be the subject of a future work. 

In this letter, we solve exactly the Fano dissipation 
problem for weak field, fulfilling the conditions for ab¬ 
sorption, Raman, Rayleigh and fluorescence spectro¬ 
scopies. Within the wideband approximation to describe 
the coupling to the continuum, we obtain explicit and 
simple expressions for all optical observables, such that 
they can be used explicitly by experimentalists, notably 
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to directly obtain the system parameters. For this, we 
solve the dynamics of the quantum system with energy 
levels in a Fano like configuration and coupled to a bath 
which induces excited states relaxation and pure dephas¬ 
ing. The system density matrix evolution is described in 
the Born-Markov approximation, by the Liouville equa¬ 
tion in Lindblad form EH- 

Although our formalism is general and can be applied 
to a variety of systems, an important motivation is the 
type of architectures encountered in light harvesting sys¬ 
tems called Gratzel or dye-sensitized solar cells EH US], 
where a molecule is adsorbed to a semiconductor sur¬ 
face. In most important physical realizations, the elec¬ 
tronic ground state of this hybrid system is an isolated 
quantum state located in the semiconductor energy gap 
while the excited states can be considered as superpo¬ 
sitions of localized molecular excited states and delocal¬ 
ized semiconductor conduction band states EQ]. Ev¬ 
idence suggests that such a model with minor modifica¬ 
tions could also account for molecules on metal nanopar¬ 
ticles na in nn- Be it for solar energy applications, 
electronics or sensors, there is strong evidence that the 
details of the interface distinguish functioning from non¬ 
functioning devices [25] , 

The energy levels of our model, along with the possi¬ 
ble transitions, are shown in Figure 1. A discrete excited 
state \e) with energy E e is coupled to a continuum of 
states \k) with energy E k . These states can be reached 
from a ground state | g) through laser excitation. A sub¬ 
manifold (typically vibrational) \v) with energy E v is in¬ 
cluded in the electronic ground state to open inelastic 
scattering channels as \v) —> \v > ). The Hamiltonian is: 


H = H 0 + H V + H f (2) 

H 0 = ' s ^E u \v){v\ +£' e |e)(e| + f dkE k \k)(k\ 
v J 

H v = j dk[V{k)\e){k\+V*(k)\k){e\] 

H f = F 

ve cos(w L f)|i/)(e| + nl e cos{uj L t)\e)(v\] 

V 

+ ^Z / dk\n uk cos(u L t)\v)(k\ + nl k cos(u L i)\k){v\\ 

V 

(3) 

where Hq is the site Hamiltonian, Hy is the coupling of 
the excited state to the continuum, for simplicity, in the 
following, we will consider that V(k) = (e\H v \k) is real. 
H F is the interaction with the incident radiation field of 
frequency ut, allowing transitions from the ground state 
to the discrete excited state v o e and to the continuum 
of states v k, i-iij = (i|/i|j) is the transition dipole mo¬ 
ment between states i and j and F is the field amplitude. 

The originality of our model consists in taking into 
account a possible population relaxation from the con¬ 
tinuum to the ground state at a rate T(fc,;/), and pure 
dephasings. When the photo-excitation involves an elec¬ 
tron, this relaxation is mainly the result of the electron- 
hole Coulombic attraction followed by the thermal relax¬ 
ation through the phonons of the environment. We phe¬ 
nomenologically capture this dissipation process using a 
superoperator of Lindblad EH form, which ensures the 
trace-preserving and complete positivity of the dynam¬ 
ical map, and solve the evolution of the density matrix 
with Liouville’s equation: 



FIG. 1: Energy levels and transitions of a Fano model with 
dissipation, including a vibrational ground state manifold, a 
discrete excited state and a continuum (bottom). A particular 
realization of a Fano system with dissipation is a molecule 
adsorbed on a metal oxide semiconductor, here two catechol 
molecules on a ( 101 ) anatase TiC >2 surface (top) 


I = m, ( 4 ) 

where C(t) = Cn{t) + L D , with HCh = —*(10 H T (t) — 
H(t) 0 1), L D = L k + L? h + Lp Ule 

L k dkT(k, v)$yA{k, v) 0 A(k, u) 

— 5 [l ® Al(fc, v)A{k, v) + A^{k, v)A{k, v) 0 l] j, (5) 

L %b = E^o r v ib{A(i/, 0) 0 A{v, 0) 

-\ [l 0 A* (u, 0 )A(u, 0) + Al (v, 0 )A{v, 0) 0 l] j (6) 

and 

^pure = -E„7 e *I|e)(e| ® \v)(v\ + k)H ® |e)(e|] 

-T.U i dk lku[\k){k\ 0 \v)(v\ + \v){v\ 0 \k)(k\] 

-/ dkj ke [\k}(k\ 0 |e)(e| + |e)(e| 0 \k){k\], (7) 

where H T denotes the transpose of H , A(i,j) = |j)(*| 
are the jump operators and T(k,u) is the population 
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relaxation rate from state | k) to \v). The superoper¬ 
ator L® h relaxes the ground states vibrational mani¬ 
fold and r y0 is the population relaxation rate within 
the electronic ground state manifold. % :l are the pure 
dephasing rates for the ij coherences. We have used 
the isomorphism LpR —> L (g> R T p, where p is the col¬ 
umn form of the matrix p through the correspondence: 
\l)(m\ O |Z) (g) \m) = || Im) [24 j. 

Non-radiative and radiative transitions from the elec¬ 
tronic excited state |e) to the ground states \u) exist. 
They have been left out for simplicity, although their in¬ 
clusion can be done straigthforwardly. Furthermore, in 
most cases of interest, the electronic coupling V induces 
an effective population relaxation rate 2imV 2 /h associ¬ 
ated to electron injection into the semiconductor band 
which is much larger than the e —> v rate corresponding 
to the direct relaxation mechanism, hence justifying this 
simplification [25] . 

The optical response (absorption and emission spec¬ 
tra) is obtained through the Fourier transform of the 
field two-times correlation function which in the far field 
is proportional to the dipole two-times correlation func¬ 
tion. Using the quantum regression theorem [25H28] . the 
emitted light differential scattering cross-section can be 
written in the steady-state of the system as [29, 301: 

dnd(Hco) ^ 

x 51 55^> 55 R>& \pra,C*a,b,rb ( ?(cu Cd^))], 

aGe,fc r£y,e,fc 

with A{0) = ft i dfl the element of solid an¬ 

gle, Ij n the incident laser intensity, 9 the polar an¬ 
gle in spherical coordinates, p the steady state density 
matrix and G(z) = (z 1 — iQ L — L) _1 the resolvent, 
which is the Laplace transform of the evolution super¬ 
operator, corresponding to the time-independent Liou- 
villian L = e l — Lt C(t)e~ l — Lt in the rotating-wave ap¬ 
proximation (RWA), which consists in removing reso¬ 
nant oscillating prefactors and discarding non resonant 
terms. The fl L matrix is a diagonal matrix with 
for excited(ground)-ground(excited) coherences, and zero 
elsewhere. In summary, in order to obtain the opti¬ 
cal emission we need to calculate the resolvent and the 
steady-state density matrix, which we do next. An ex¬ 
plicit expression for G(z) is obtained by separating L into 
its diagonal and non-diagonal parts. While the diagonal 
contribution leads to a trivial calculation of its resolvent, 
the non-diagonal contribution is worked out through a 
Dyson equation to all orders in the field-free coupling 
and to second order in the field interaction. As for the 
steady-state density matrix p , it appears as the kernel of 
the Liouvillian within the RWA, that is (— iQ L — L)p = 0, 
which can in turn also be expanded to second order in the 
field, with the help of the previously calculated resolvent 

G(z) [25]. 


In principle, the above method allows to calculate the 
optical response for arbitrary couplings and density of 
states. Nevertheless, for many materials, the wideband 
approximation, which considers fc-independent couplings 
and a constant density of states, can faithfully reproduce 
experimental measurements. This was the same approx¬ 
imation which allowed a closed analytic form in the orig¬ 
inal Fano model [3] and which has been used in all previ¬ 
ous partial attempts to solve the Fano-dissipative prob¬ 
lem dflSI. We introduce the notation T cl/ = T(k,v) 
and p, vc = p u k for the /c-independent parameters. 

A straightforward but tedious calculation [25] yields 
all the terms needed to calculate the absorption spec¬ 
trum and the emission differential cross-section (Eq ([ 8 ])). 
The absorption and emission profiles, which carry the u>l 
dependence, can be expressed in a compact and simple 
form in terms of the following function: 


f{e,q,r],a) 


= a 


(g + e ) 2 

e 2 + 1 


+ V 


g + 1 
e 2 + 1 


(9) 


which is a linear combination of a Fano profile (first term) 
and a Lorentzian (second term); q being the Fano asym¬ 
metry parameter and e the normalized incident laser en¬ 
ergy. a (which can only take the values 0 and 1) and 
p are weighting coefficients for the Fano and Lorentzian 
components, respectively. 

The sum of excited state populations ( p ee + / dkpkk) 
can be written as: 


Aexcited = £ ab 7( e ,g,77,a) (10) 

where £? abs is a constant [25], e and q are the stan¬ 
dard Fano parameters (e = (E e — 7 = mrV 2 , 

q = poe/{mrVpoc)), a = 1, and tj = r cJ ,)/(47mV 2 ). 

We can obtain the extinction coefficient ft measured in 
optical absorption measurements using arguments of de¬ 
tailed balance 1 25] : 


ft = n7T ^° c hu L f(e, q, rj = 0, a = 1) (11) 

ceo n 

where c is the speed of light and eo the vacuum permittiv¬ 
ity. The differential cross sections for the three processes 
Rayleigh, Raman and fluorescence can all be expressed 
referring to the same functional form: 


dV 

d£ld(Hu}) 


A(9)BlR(u J ,u J iA i )f(e i ,q i ,r ] i ,a i ) 


( 12 ) 


where i = Rayleigh, Raman or fluorescence. The ex¬ 
pression for emission has 3 parts: A(9)Bl is a prefac¬ 
tor [25], R(uj,ujo, A) = ^[(w — wo) 2 + (A) 2 ] -1 is a normal¬ 
ized Lorentzian of central frequency ujq and half-width A 
which gives the emission lineshape (for a fixed incident 
laser frequency), and f(e,q,r),a) which gives the profile 
- the signal integrated over the emitted frequency - and 
carries all the dependence of the incident laser frequency 
WL- 
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We first discuss the characteristic function / of the 
population given by equation 10 The result is the com¬ 
bination of a standard Fano profile, and a Lorentzian. 
The signature of the relaxation process is embodied in 
the parameter 77 . Its physical meaning is the ratio of the 
relaxation rate to the injection rate into the continuum. 
When 77 1 the relaxation quenches the population of 

the continuum and therefore suppresses the Fano inter¬ 
ference signal giving rise to a pure Lorentzian lineshape. 
When i] <C 1, we recover the original Fano profile. Fig¬ 
ure [ 2 ] sh ows normalized population profiles as given by 
Eq. (|10[), for several values of q and 77 . As 77 —> 0, the 
curve approaches the standard Fano profile while 77 = 1 
shows the Fano with a Lorentzian contribution for differ¬ 
ent values of the q parameter. 



FIG. 2: Characteristic function for absorption or Raman 
emission for 77 = 0 (solid black) and 77 = 1 (red dashed) for 
q = 0 , 1,100 (intensities are normalized). 

For the emission processes, we separate the full opti¬ 
cal response given by Eq. ([8]) into several components 
according to their central emission frequency Uq of the 
lineshape functions R(w,Wg,A l ). The Rayleigh process 
corresponds to uJq ay = u>l and the Raman process to 
w Ram _ Ul _ Wvibj (IicUvib = Ei — Eq). Contrary to 
the Rayleigh and Raman scattering which are coherent 
processes, fluorescence corresponding to a radiative tran¬ 
sition from excited population is incoherent. This sta¬ 
tionary population can only exist if in addition to pop¬ 
ulation relaxation, pure dephasing processes are taken 
into account (777 7 ^ 0 in Eq. 0 ). The population of 
the discrete excited state leads to an emission at E e /h 
and at E e /h — w v ib central frequencies. The contin¬ 
uum states are populated at the laser frequency uil and 
therefore give a fluorescence emission with central fre¬ 
quencies at Wi and wl - w v ib- We note that the fluo¬ 
rescence from the continuum states occurs at the same 
central frequencies as the Rayleigh and Raman scatter¬ 
ing, but can be distinguished by their respective line- 
shape widths A*. For Rayleigh, the width comes from 
the laser linewidth, A Ray = <5, and for the Raman the 
width is given by the inverse lifetime of the vibrationally 
excited state v = 1, A Ram = r v ;b/2. The width of the 
fluorescence lineshape is dominated by the excited state 
population lifetime: for the emission at u>l the width is 
^fluor,cO _ Y cu and for the emission at uil — UMb the 
width is A fluor ’ cl = T cv + T v ib/2. All the parame- 
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FIG. 3: Comparison between different Fano models. Top: fits 
of the profile presented in this paper Eq. with a) q = 4 and 
b) q = 1 , with r) = 1 by a standard (blue dashed) or a shifted 
Fano profiles (solid green). Bottom: relative differences be¬ 
tween the extracted f/ e ff and the actual q value for c) q = 4 
and d) q = 1 for the standard (blue dashed) or the shifted 
Fano profiles (solid green). 


ters for each emission component given by Eq. (12 1 are 
collected in Table HI 

The parameters of the model can be extracted from 
spectroscopic measurements. By fitting the profiles, we 
can get 7 , q , R abs , £? Ram , _B Ray , F (experimental parame¬ 
ter), 77 Ray , and 77 Ram (Tableland S.M. [25]). From these, 
there are various ways to obtain the model parameters, 
for instance: 


^ 5 Ray F 2 ^ 

(£ »ab s) 2 2h' ^ V ~\J n 
I £Ray / B Ray 7 

= V Moe = V 

g^Ram^TjabsTjRam^ 


(13) 


_ ' ^Ram 

c = \/ 7r ^ a b s » A*le = 


77 ffR a y F 2 


We now discuss how important is to include the 
Lorentzian term arising from dissipation in addition to 
the standard Fano profile. In Fig. [3] we present a profile 
with our model Eq. (10 1 with a given set of parame¬ 


ters and fit it both with a standard Fano model /i(e; q e ff) 
(Eq- 0 ) and with a shifted Fano model q e g)+D], 
The shifted Fano fitting is usually performed in optical 
measurements and corresponds to a constant background 
substruction. Below, we calculate the relative error be¬ 
tween the extracted value of the Fano asymmetry param¬ 
eter, which we call q e g, and the original q as a function of 
the Lorentzian weight factor 77 . We show that even when 
the fitting is very good, the extracted parameters can be 
off by a factor of two for 77 = 1. This definitely shows the 
importance of including properly the dissipation in the 
Fano model. 

In conclusion, we expect that the equations will mo¬ 
tivate the experimentalists to extract system parameters 
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Lineshape 

Profile 

Model: 

Center 

frequency 

Width 

Laser Asymmetry 
frequency parameter 

Fano Lorentzian 

weight a weight 77 

Standard: 



e q 

1 

Excited state: 

populations 



e q 

1 77 

Rayleigh: 

Raman: 

ojl 

OJl — Wvib 

5 

Fvib /2 

e q 

e Q 

1 -^V?7 

Oc 

1 -^7? 

Fluorescence: 

discrete 

(E e - E 0 )/h 

(E e -E 1 )/h 

n7rV 2 3 4 5 6 7 8 /h + 7e0 
nnV 2 /h + 7 e i + r vib /2 

e(mrV 2 ) q(mrV 2 ) 

mrV 2 +^7e0 mrV 2 +^7e0 
e(mrV 2 ) q(mrV 2 ) 

mrV 2 +ft'Y<=n nirV 2 |^ pn 

n ( n Z v ‘ 1 ) 2 n 

k mrV 2 -\-h , y e Q ' • 

n ( mrV 2 \2 

' mrV^-\-h'y P o ' ' 

Fluorescence : 

continuum 

OJL 

OJL ~ Wvib 

5^1/ Tczy + 2 7e0 

T„. r Cl/ + 7e0 + 7el + Tvib/2 

e(mrV 2 ) q(mrV 2 ) 

1 TeO 1 

mrV 2 +/i7 e o 7i7rV^+fi^y e o 
e(mrV 2 ) q(mrV 2 ) 

(mrV' 2 /h+ 7 e0 y 2 q 2 +1 

1 ^eO 1 

7i7rV^+/l7 e o mrV 2 -\-fr'YeO 

(mrV 2 /h-\-'Y e Q) 2 q 2 + l 


TABLE I: Parameters for all optical processes depicted in Eqs. (10 1 (12 1 , 
e = (Swl - E e )/(nnV 2 ), q = p 0e /(nnV p 0c ) and 77 = J2„ T' cv /(innV T ). 


The parameters are expressed as a function of 


which were not accessible previously, from routine op¬ 
tical spectroscopies, characterizing with increasing pre¬ 
cision the discrete-continuum interface both relevant in 
devices and interesting from a fundamental standpoint. 
Furthermore, our model serves as a stepping stone for 
further theoretical studies: optical response beyond the 
wide-band approximation (near band edges for example), 
strong field effects relevant for plasmonics, generalization 
to time-dependent laser pulse sequences relevant to non¬ 
linear 2D spectroscopy and application to real systems 
with the help of computational tools in order to obtain 
ab-initio the parameters of the model. 
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